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Abstract 

The trace of a function of a Schrodinger operator minus the same for the 
Laplacian can be expressed in terms of the determinant of its scattering ma- 
trix. The naive formula for this determinant is divergent. Using a dispersion 
relation, we find another expression for it which is convergent, but needs one 
piece of information beyond the scattering matrix. Except for this 'anomaly', 
we can express the Casimir energy of a compact body in terms of its optical 
scattering matrix, without assuming any rotational symmetry for its shape. 
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1 Introduction 



The kinetic theory of gases in classical physics predicts that every surface 
in a gas is bombarded by molecules. The recoil of these molecules exert 
a pressure on the surface. A spectacular demonstration of this prediction 
would be to evacuate the air inside an aluminium can: there will be a net 
inward pressure that will cause the can to collapse. 

In the quantum theory of fields, there is an analogous pressure on every 
surface that can scatter light [1, 2]. Even in the vacuum there are 'virtual 
photons' which describe the quantum fluctuations of the electromagnetic 
field in its ground state. These virtual photons are scattered by any medium 
that can interact with light; this scattering exerts a force on the medium. 
Although quite small in magnitude, it has been measured experimentally [3]. 

In the original calculation of the Casimir force only simple shapes such 
as a sphere, or a rectangular slab were studied. Moreover, the medium was 
assumed to have ideal properties, such as perfect conductivity. Since that 
time, new methods have been developed, which allow the calculation of the 
Casimir energy for more general shapes. Of particular interest to us are the 
spectral methods[4, 5, 6, 7, 8] . (Methods that allow realistic computation 
of Casimir energy of micro-electromechanical or MEM devices have been 
developed recently [9]; but we do not use them in this paper.) The basic 
idea behind these methods is a relationship (often called Levinsion's theorem 
or Krein's trace formula) between the optical scattering matrix S and the 
density of states: 
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The Casimir energy is the sum over frequencies weighted by p(k). 

The physical argument that the Casimir force is due to the reflection of 
virtual photons suggests that it should be possible to express it entirely in 
terms of the probability of reflection. Unless the momentum of the photon 
changes during scattering, it should not contribute to the Casimir force. 

In this paper we will show that the dispersion relations of scattering the- 
ory allow an answer for the density of states in terms of the reflection proba- 
bility alone, in the one dimensional case. In the three dimensional case, there 
is a potential logarithmic divergence in tr log S . We show how to remove 
this by using a dispersion relation, without assuming rotational symmetry 
or using perturbation theory as in previous treatments [10]. An outcome is 
an 'anomaly' term in the density of states proportional to the integral of the 
potential. Contrary to the physical picture in terms of scattering of virtual 
photons, there is this one contribution to the Casimir effect that cannot be 
expressed in terms of the scattering matrix alone. 

Much of our detailed analysis will be only carried out for the technically 
simpler case of a scalar field. As in traditional optical scattering theory, 
this scalar model gives a reasonable picture of the essential phenomenon. 
It is straightforward although technically involved to extend our analysis 
to include polarization effects; i.e., scattering of vectorial fields. Moreover, 
we will ignore the effects of absorption of light: the effect of absorption on 
virtual photons needs a deeper analysis than we can provide at the moment. 
Also, we will consider only optical media that are time independent. For a 
technical reason, we will further assume that the scatterer is parity invariant; 
although it does not have to be symmetric otherwise. 
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We aim to give a more or less self-contained description, taking off from 
the discussion of scattering in standard textbooks. The erudite reader could 
skip ahead to sections 3.1,3.2 and 5. No pretense at mathematical rigor is 
made. However, there will be an occasion for us to be careful about the 
definition of an infinite determinant (of the scattering operator) to avoid a 
logarithmic divergence. 

2 Scattering of Light 

The propagation of light of(Ref. [11]) wave-number k in a dielectric medium 
is described by the equations 



Here, e(x, k) is the dielectric 'constant' which may in fact depend on 
position and on the wave-number k. (In other words, e(x, k) is the square of 
the refractive index.) E and H are complex vector valued functions describing 
the amplitude of the electromagnetic wave. 

Eliminating H, we get 



VxE = ikH, V x H = -ike(x, k)E. 



(2) 



V x V x E(x, k) = k 2 e(x, k)E. 



(3) 



In other words, 



V x V x E(x, k) + V(x)E = k 2 E 



(4) 



where the 'effective potential' [12] is 




(5) 
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There is an analogy between the above equation describing the classical 
scattering of light and the quantum mechanical scattering of a particle by a 
potential. The Schrodinger equation (in natural units Ti = 2m = 1) is 



The essential difference is that the vector Laplacian V x V x E is replaced 
by the scalar Laplacian V 2 . The potential V(x) is related to the refractive 
index as in the formula above. If the refractive index is less than one, the 
'potential' is positive; this is the case for large enough k. 

Thus it will be useful to use the language of quantum mechanical scat- 
tering theory, although of course the scattering we are talking about is just 
the classical scattering of light waves. But some differences must be kept in 
mind. For example, the 'potential' V(x, k) will depend on wavenumber, a 
fact that has no simple meaning in quantum mechanics. The wave-number 
dependence of the scattering matrix will be affected by this dispersion of the 
light waves. Our final formulas only involve the scattering matrix and thus 
take into account of this effect. But we will often not explicitly display the 
dependence of the potential on the wavenumber. 



We will start our discussion in the simplest possible case, the scalar theory in 
one space and one time. We will later generalize the ideas to three dimensions 
and to vector waves, but it is pedagogically useful to start with the simplest 
case. 



V 2 ifj + V(x)i/> = k 2 ij. 



(6) 
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Waves in 1 + 1 Dimensions 
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Often we are interested in deriving a formula for the sum over a function 
of the eigenvalues of a hamiltonian operator H. This can be thought of 
as the trace of an operator tr <f)(y/H). For example, the Casimir energy of 
a massless scalar field in the presence of an external potential is tr y/H. 
Another example is the thermodynamic partition function, tr e~^ H . More 
precisely we will be interested in the difference between this trace and its 
value for the free hamiltonian. 

With this in mind, let 

We can think of H be as the Schrodinger operator of a one dimensional 
quantum mechanical system. The potential V(x) is positive, smooth and 
vanishes faster than any power of x at infinity. Thus, H > 0. 

Let <p(k) be a smooth function on [0, oo) which vanishes at infinity faster 
than k~ x . We are going to derive a formula for the quantity 

tr <j>(y/H) = tr [</>(y/H) - (f>(y/H )]. (8) 

We assume, for now, that <j>{k) vanishes for large k : postponing the study 
of ultraviolet divergences. We might expect that there is an integral repre- 
sentation 

POO 

tr 4>UH) = / <f>(k)p(k)dk. (9) 
Jo 

We will get a formula for the 'spectral density' p(k) in terms of the scattering 
matrix. We can get this directly from the formula for the resolvent in terms 
of the Jost functions[5, 6], but we will give a more physical but less rigorous 
argument. 
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Instead, we will look at the continuous spectrum as the limit of an eigen- 
value problem. Imagine we have enclosed the whole system in a box of size 
2L: we require that the eigenfunctions of H and H vanish at x — ±L. 
Then we have a discrete number of allowed values of momentum, k n (L). For 
example when V(x) = 0, they are |g. Thus 

tr ^H) = J2 

n 

This sum will become an integral in the limit as L — > oo. Even when L 
is kept finite, we can assume that it is much larger than the range of the 
potential V(x), since we are only interested in the limit L — > oo. 
At infinity, the solutions of the Schrodinger equation [13] 

-u" + V(x)u{x) = k 2 u(x) (11) 

tend to a sum of plane waves: 

u(x) Ae ikx + Be~ ikx for x -oo 

-> Ce lkx + De~ ikx for x -)■ oo. (12) 

Thus A, D are the amplitudes of the incoming waves and C, B those of the 
outgoing waves. Conservation of probability gives 

\A\ 2 + \D\ 2 = \C\ 2 + \B\ 2 . (13) 

Since the Schrodinger equation is a second order differential equation, there 
are two constants of integration in its general solution. Let us choose them 
to be A and D, the amplitudes of the incoming waves. Then the amplitudes 



<t>{h n {L)) - (f>( 



7rn, 
2L' 



(10) 
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of the outgoing waves are determined once we solve the differential equation: 

The matrix S(k) is unitary from the conservation of probability: it is the 
scattering matrix associated to the potential V. 

Often we describe the scattering in terms of reflection and transmission 
coefficients. They are defined by the special choice A — 1, D — correspond- 
ing to waves incident from the left: 

u(x) ->• e ikx + R{k)e~ lkx for x -> -oo 

-)• T(A;)e ifcc for x ->■ oo. (15) 

Taking complex conjugate we have another solution, 

u*(a;) ->• e- ifee + ir(*;)e ifca! fora;->-oo 

->■ T*(k)e- ikx for x ->■ oo. (16) 

Taking linear combinations, gives the solution describing waves incident from 
the right: 

u*(x) - R*(k)u(x) 1 - ^(A;)! 2 _., 

— ' ,,\ ->■ — ' e * fe:r for x ->■ -oo 

T*(k) T*{k) 

^ _ R*(k)T(k) ikx e _ lkx for 

T*(jfc) V ; 

Thus (recalling that \R(k)\ 2 + |T(A;)| 2 = 1), 

(T{k) -R*(k)f§r\ 

S(k) = 18 

We might check explicitly that this matrix is unitary. 

8 



Being unitary, there is a basis in which S(k) is diagonal, with eigenvalues e 2M?1 > 
The real valued functions 771,2 (k) are the 'phase shifts' of the scattering prob- 
lems. 

Now let us return to the eigenvalue problem 

-u" + V(x)u(x) = k 2 u(x), u(-L) = u(L) = (19) 

Since L is large compared to the range of V(x), we can use the asymptotic 
forms above in the boundary conditions: 

Ae -ikL + Be ikL _ — Ce ikL + De~ ikL . (20) 
Solving for B and C, 

B = -Ae~ 2ikL , C = -De' 2ikL . (21) 
Thus the momenta k are determined by the equation 

S{k) {t) = ~ e ~™ L [a) ' (22) 

So far we are in arbitrary basis in the two dimensional space of solutions. 
Now let us choose this basis to the one that digaonalizes S(k): 

e 2im(k) A = _ e -2ikL D ^ e 2i V2 (k) D = _ e -2ikL A ^ 

Eliminating D, 

e 2i[r)i(k)+r) 2 (k)+2kL} _ ^ ^4) 

This is the transcendental equation for the allowed wavenumbers, 

t = g- ^+^*) ,fc,n= ■■-2,-1,0,1,2.- (25) 
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Since L is very large in the limit we are interested in, the second term is 
a small correction. It is enough to solve the equation approximately by 
iterating it once: 

(26) 



We can now reexpress the sum over phase shifts in terms of the Scattering 
matrix: 

2i[r )l (k)+r ]2 (k)} = logdetS(£;). (27) 
For later use we also note that, 



detS(fc) = T 2 (k) + 



2 , t , , \R(k)\ 2 T(k) _ T(k) 



so that 



Thus, 



T*(K) T*(k) 
logdetS(A;) = 2iaigT(k). 



nn_ logdet^(fg) 
ny ' 2L 4iL 



(28) 



(29) 



(30) 



Now we are ready to evaluate the sum over momenta in the limit as 
L — > oo. 



(t>{k n {L)) - 0( 



7m, 
2L' 



^ nn logdetS^ 
-? 0( 2Z } - 



Trn\ 
2L> 



AiL 



(31) 



This becomes an integral as L — > oo: 



(32) 



2m d/c 

Similar trace formulas have been derived by many other methods[5, 6, 10]. 
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If the potential is independent of frequency, the integral will converge if 
4>{k) falls off faster than k^ 1 for large k and remains finite for small k. (This 
is because T(k) ~ 1 + ^ / V(x)dx for large from the Born approximation; 
more on this later.) If <p{k) does not fall off at large k, the integral can still 
converge if V(x,k) tends to zero for large \k\, so that S(k) tends to unity 
faster also. 

Even if the potential depends on the wave-number, it will still determine 
a unitary scattering matrix S(k) for each k. We can still ask for the sum 
over the allowed values of wave-number for finite L as before; it is determined 
by the same formula in terms of the scattering matrix. Indeed, the whole 
argument goes through without any change. 

3.1 Dispersion Relation 

We will now express the density of states in terms of the reflection coefficient. 
The force due to the scattering of a virtual particle ought to be dependent 
on the probability of reflection: if the particle is transmitted, its momentum 
does not change and hence it should exert no force. Thus it is physically 
unsatisfactory that we have a formula in terms of the argument of the trans- 
mission amplitude rather than the magnitude of the reflection amplitude. 

But we now remember the basic fact that T(k) is analytic in the upper 
half plane. There are no poles since there are no bound states. Moreover, 
T(k) never has zeros in the upper half plane. Hence logT(fc) is also analytic 
in the upper half plane. There is then a dispersion relation [13] between its 
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real and imaginary parts (Hilbert transform) : 



But 



im bgrw = / "^ "W . (33) 

7r J k — k 



Re logT(fc) = ilog|T(A;)| 2 = 1 log[l - |i?(£)| 2 ]. (34) 



Thus 

argdet5(Jb) = 2argT(Jb)=P / !^df)3 g, (35 ) 



3.2 The Trace In Terms of the Reflection Coefficient 

Combining the above results we get 

&<MH) = -/°°^^argdet5(fc) (36) 
Jo Zir ak 

p r*r d^ ^w logji-i^)! 2 ] , 37) 

Jo 27ri-oo 7T c?fe k-k' { ' 

Assuming parity invariance |i?(fc)| = \R(— k)\, 

We can expand the log in a power series to get a 'multiple reflection expan- 
sion' for this quantity. 

When (f){k) = k this can be simplified to 

1 r dkdk> klog[l-\R(k'^]-k'log[l-\R(k)\^] 
tlVH ~ 2 Jo W^m (39) 

with a non-singular integrand. 
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4 Waves in 3 + 1 Dimensions 



Consider the operator H = — V 2 + V(x), on L 2 (R 3 ) where the 'potential' 
V(x) is a real function on R 3 that vanishes at infinity. (The potential also can 
depend on the wave-number: but again we don't display this explicitly.) We 
will also assume for simplicity that that H has positive continuous spectrum 
; i.e., that it has no 'bound states'. 

We will again be interested in the the difference 

tr(0(v/#)) = tr <j>{y/H) - tr 0(^o) (40) 

of the trace from that with the free hamiltonian: 

H = -V 2 . (41) 

The function (f>(k) will, to begin with, be assumed to be smooth and vanish at 
infinity faster than A; -3 : this is necessary to avoid 'ultra-violet' divergences. 
Later on we will allow for (f>(k) = k and check that the integrals converge 
when the wave number dependence of the potential is take into account. 
Again we will derive a formula that involves the scattering matrix. 

Since the potential vanishes at infinity, there is a solution to the equation 

[-V 2 + V(x)\tj)(x) = k 2 i)(x) (42) 

that tend to a sum of a plane wave and an outgoing spherical wave [13]: 

p ikr 

Mr, n') ->• e ikrn - n ' + f(k, n, n') , as r ->■ oo (43) 

r 

Here, n is the direction of the incoming plane wave and n' the direction along 
which we let the argument of the wavefunction go to infinity. The function 
f(k, n, n') is the scattering amplitude. 
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Since there are no bound states, the general solution can be written as a 
superposition of such scattering solutions, [13] weighted by a function F(n) 
of the direction of incidence. Such a general solution will have the asymptotic 
behavior at spatial infinity, 

Mr, n') ~ / F(n)e ikr ™'dn n + ^ / F(n)f(k, n, n')^„ (44) 

As r — > oo the first integral can be evaluated by the method of steepest 
descents, to get 

p—ikr ikr 

2*iF(-n')—-2mF(n') — . (45) 

Thus we can write 

2tp? r p~ ikr p ikr ~ ~) 

M^ n') ~ -Pi ! F (- n ') 1 1 + 2^/(fc)F](n') (46) 

k I r r J 

where /(&;) is the operator on L 2 (S 2 ): 

/(fc)F(n') = y /(fc, n, n')F(n)^ n . (47) 

The quantity 

S(ife) = 1 + 2ikf(k) (48) 

is the scattering operator. Conservation of probability requires it to be uni- 
tary. 

Being unitary this operator can be diagonalized on L 2 (S 2 ). If H is spher- 
ically symmetric, the eigenf unctions of S are the spherical harmonics: 

S(k)Y lm (n) = e 2 ^Y lm (n). (49) 

The eigenvalues are determined by the 'phase-shifts' rji(k) in each angular 
momentum sector. Even if H is not spherically symmetric, there will be 
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some spectrum of 'phase shifts' r) a (k): 

S(k) Xa (k,n) = e 2 ^Xa(k,n). (50) 

(We will assume for now that this spectrum of S is discrete; i.e., that the 
label a takes values in a countable set.) 

Let us first imagine that our whole system is enclosed in a spherical box 
of large radius R with the wavefunction required to vanish on this surface. 
(The exact shape doesn't matter in the limit R — > oo, which is all we are 
interested in, so we might as well assume it is spherical.) For finite R there 
is a discrete set of allowed values for k, say k na (R). They are fixed by the 
condition 

V>CR,n') = 0. (51) 

When R is large, this becomes, 

S(k)F(n) = e- 2lkR F(-n). (52) 

By squaring, 

S(k) 2 F(n) = e- 4ikR F(n) (53) 

Not every solution of (53) is a solution of (52). In the limit of large R, we 
can take account of this over-counting just dividing the density of states by 
two 1 . 

^or example, in the spherically symmetric case each solution is degenerate, with de- 
generacy labelled by the angular momentum quantum numbers. In order to satisfy (52), 
we have the additional condition that n is odd for odd angular momentum I and even for 
even I. This is not needed for (53) 
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The solutions of (53) are the eigenfunctions of the scattering matrix 
introduced above; the momenta are fixed by 

Ar) a (k) + AkR = 2im, ra = 0,l,2,--- (54) 

Thus, the solutions are labelled by n, a and will depend on R: k na (R). In the 
case of the free particle, the phase shifts vanish and the allowed momentum 
values are 

As R becomes large the solutions will differ from this by a small correction: 

M«) = ^ + ^>. (56) 

We get 

7TT? 

q na {R) = -Va(^). (57) 
Now consider the difference of traces, (remembering to divide by two to 
avoid the above mentioned over-counting) 

tr 4>WH) = \Y\4>{k na {R)) - ■ (58) 

As R — > oo we get 

^WH) = -^/(^Q - -Ef * (59) 

In other words, 

tr 0(Vif) = | rffc^ tr log5(fc). (60) 

We are thinking of S(k) as an operator on L 2 (S 2 ), so the trace on the r.h.s. 
is the average over angles. We will now turn to a more precise definition of 
this trace: naively it is logarithmically divergent. 
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4.1 A toy model: the Gamma function 



The reader familiar with the theory of modified determinants as in equation 
(67) below could skip this subsection. 

Recall that the Gamma function has poles at negative integers and zero; 
its reciprocal is an entire function with zeroes at these points. So we might 
hope for a product formula 

1 



T(z) 



n 

n=l 



n 



(61) 



Alas, this product is divergent: the product l\ n [1 + \i] converges when the 
sum J2 n \ converges. In the our case, J2 n \ diverges logarithmically. But, 
the sum of the squares converges: J2 n ^ < 00 • 

This suggests a fix to this divergence problem. The function e~ 2 (l +z) has 
the same zero as 1 + z but tends to one faster as \z\ — > 0. 

Indeed, |e~ z (l + z) — 1] < C|z| 2 .So the product 



n 

n=l 



e n 



n. 



(62) 



converges. So we can separate out a divergent part of (66) as follows: 

OO 1 



1 °° f _* 

log — — = log z + lo S 1 e n 









je " 




} 




n_ 





(63) 



n=l 



The divergence has been isolated to the last term. How to give a meaning to 
that last divergent sum? We might guess that the correct formula is 
1 



log 



log z + J2 1°§ \ e 



z 

1 + - 

n. 



+ az 



(64) 



r(^) 

for some constant a; and determine it by comparison with the value of the 
logarithmic derivative of T(z) at some point. This is a kind of renormalization 
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of the logarithmic divergence. Indeed, we know that 



d 



iogr(z) 



1 



7 + 0(z) 



(65) 



dz 



z 



where 7 is the Euler constant. Thus we conclude that = 7: 



1 



00 




ze fZ J] 



n=l 



e « 1 H — 
n. 



(66) 



This argument would not pass muster with a modern analyst. But it is 
precisely such heuristic arguments that led to the rigorous modern theory of 
analytic functions. Quantum Field Theory is still in the stage of development 
that complex function theory was in the mid nineteenth century: we need 
to work with heuristic, physically motivated arguments which point the way 
to the truth. These then should be turned into theorems later, as stronger 
constructive methods become available. 

4.2 Infinite Determinants 

The determinant of an operator 1 + X is well-defined when X is trace-class: 
that is, when the sum of the absolute values of its characteristic values is 
convergent. But the scattering matrices of interest to us are not this type: 
S — 1 is not trace-class. However, we will see soon that it has the next 
best property: the sum of the absolute squares of the characteristic values 
of S — 1 converges, because it is proportional to the total scattering cross- 
section. WhenX is Hilbert-Schmidt (i.e., tr X*X < 00), the modified 
determinant [14, 15] 



deti [l + X} = det e~ x [1 + X] 



(67) 
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is meaningful. (The point is that the analytic function e~ 2 (l + z) — 1 is 
bounded by c\z\ 2 for some constant c. Hence e~ x [1 + X] — 1 is trace-class 
if X is Hilbert-Schmidt. ) Now, we can write the logarithm of the original 
determinant as 

log det [1 + X] = log deti [1 + X] + tr X (68) 

which separates out the divergent term tr X. 
Lemma S(k) — 1 is Hilbert-Schmidt 

Proof: In terms if the scattering amplitude, 
tr (S(k) - l)t(S(fc) - 1) = 4k 2 tr pf = 4k* J \f(k, n, n')| 2 ^^ (69) 
But the total cross-section for a beam incident along the direction n is 

a(k,n) = J \f(k,n,n')\ 2 dn n ,. (70) 

Thus 

tr (S(k) - l)\S(k) - 1) = -k 2 a(k) (71) 

7T 

where 

a(k) = ^- J dQ n a(k, n) (72) 

is the average total cross-section, which is finite as was promised. 

Then (1 + 2ikf)e~ 2lkf — 1 is a trace-class operator and the modified de- 
terminant detiS(k) defined by [14] 

deti[S(£;)] = det[(l + 2ikf) e - 2ikf ] (73) 
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exists. Moreover we have the bound, 

| logdet^fc)]! < \ tr [S(k) - l]t[5(fc) - 1] = ^k 2 a(k) (74) 

which will be useful later. 
We can write 

log det S(k) = log deti5(fc) + 2iA; tr /(jfc) . (75) 

The real parts of the terms on the l.h.s must cancel, since the scattering 
matrix is unitary. So we might as well write, 

log det S(k) = i Im logdet^A;) + 2ik Re tr f(k). (76) 

5 A Dispersion Relation 

The last term, which is potentially divergent, can be given a meaning using 
the dispersion relation which relates it to the scattering cross-section: (See 
[13]) 

1 r 1 r°° 2k' 2 rr(k') 

Re tr /(*) = -- / V {X , t) A + / -j^O/. (77) 

(Here, V stands for the principal part of the integral.) The first term is just 
the Born approximation of the forward scattering amplitude. We will see in 
a minute that 

a(k) ~ k~ 2 , as k ->■ oo. (78) 

(This is for the case that the potential V(x) is independent of k. It falls off 
even faster if the potential itself vanishes for large k.) Hence this integral is 
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convergent. Thus we have the formula we seek: 



~ k f k r°° k /2 (r(k f } 

argdet S(k) = arg detiS(fc) - - / V(x, k)d?x + —V / ;/9 \t dk'. (79) 



7T 



7r 2 jo A;' 2 — k 2 



As in (37) we can now get trace of a function of the hamiltonian: 

dk d(f){k) k 



fC 

tT<j>{y/H) = / 
J 



/ 

Jo 



2% dk 7r 
*> dk d<p(k) k 



J V(x, k)d} 



x 



2vr dk 7T 2 Jo k' 2 -k 2 
dk d(f){k) 



o 2-7T dk 



-arg det 1( S(fc) 



(80) 

(81) 

(82) 
(83) 



The last term, which is the most complicated to calculate, can be bounded 
(74) by the total cross-section : 



f 

Jo 



dk d4>{k) 
2tt dk 



arg det 1 S(k) 



< 



f 

Jo 



°° dk d<P(k) k 2 a(k) 



2vr dk 2vr 



(84) 



5.1 Application to Casimir Energy 



An application of the above formula is to the Casimir effect, 4>{k) = k. To 
avoid divergences, we assume that 



for large k. 



k 2 V(x,k) ->■ 



tiy/H = r^h. r v(x,k)d 3 x 

Jo 2% 7T J 

r™ dkdti_ ktt_ k'a(k') - ka(k) 

'Jo 2^2^ ~V k' 2 - k 2 

f°° dk , . 
- / —wgdet^ik) 
Jo 2ir 



(85) 



(86) 
(87) 
(88) 
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The first two terms should give a good approximation in many cases. Again, 
the last term can be bounded by the cross-section, 



The computation of the scattering matrix from the optical potential is 
a separate problem: it can be done numerically. Or we could measure the 
scattering matrix experimentally and use it directly to calculate the Casimir 
energy. 

Despite the physical intuition that the Casimir force is due to reflection of 
virtual photons, we do not get an answer entirely in terms of the scattering 
matrix: the first term involves the potential itself. This is similar to the 
way that the product formula for the Gamma function does not just involve 
the location of its poles: our term involving the potential is the analogue 
of the term involving the Euler constant in the formula for logT(z). Such 
'anomalies' occur elsewhere in Quantum Field Theory[16] as well. 
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